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Abstract 
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' The measurement of the polarized gluon distribution function AG/G using current- 

[~^. . target correlations in the Breit frame of deep inelastic scattering is proposed. The 

approach is illustrated using a Monte Carlo simulation of polarized ep-collisions for 
, HERA energies. 
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1 Introduction 



A direct determination of the polarized gluon distribution is of importance for under- 
standing the spin properties of the nucleon. Recent experimental measurements of spin- 
dependent structure functions [|]-|| have shown that the valence quarks account for only 
a small fraction of the nucleon spin. One of the possible explanations for this observation, 
known as "spin crisis" , is to assume a large contribution from the gluon spin. There exist 
many theoretical models which explain this phenomenon, but there is no experimental 
evidence to favor one of them. 

A direct way to solve this puzzle is to measure the gluon density in polarized lepton- 
nucleon scattering. One of the suggested methods is based on the detection of the corre- 
lated high-pt hadron pairs in polarized deep inelastic scattering (DIS) j|. This measure- 
ment has already been performed by the HERMES Collaboration ||. Another possibility, 
which is planned to be studied at the COMPASS experiment, is to analyze events with 
open/closed charm ||. At the HERA ep collider, after possible upgrade to polarized 
beams, the polarized gluon density can also be measured from dijet events |0]. 

In this paper a possibility to measure AG/G from the multiplicity correlations in 
neutral-current DIS is discussed. The proposed method is based on the Breit frame ||. For 

1 On leave from the Institute of Physics, AS of Belarus, Skaryna av.70, Minsk 220072, Belarus. 
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Figure 1: A schematic representation of the Breit frame for the quark-parton model. 



this frame, in the quark-parton model (QPM), the incident quark carries Q/2 momentum^ 
in the positive Z-direction and the outgoing struck quark carries the same momentum in 
the negative Z-direction. The phase space of the event can be divided into two regions 
(see Fig. [l]). All particles with negative p^ rcit components of momenta form the current 
region. In the QPM, all these particles are produced from hadronization of the struck 
quark. Particles with positive p® rmt are assigned to the target region, which is associated 
with the proton remnants. 

For the QPM in the Breit frame, the phase-space separation between the particles 
from the struck quark and the proton remnants is maximal. Therefore, it is expected 
that there are no correlations between the current-region particles and particles in the 
target region of the Breit frame |J. Such a separation does not exist when the leading- 
order QCD processes, known as Boson-Gluon Fusion (BGF) and QCD-Compton (QCDC) 
scattering, are involved JTO] . The kinematics of these processes lead to current-target anti- 
correlations predicted in |J] and experimentally measured by the ZEUS Collaboration [ 11 1 . 
The magnitude of these correlations at small x is mainly determined by the BGF processes. 



2 Current-Target Correlations. Unpolarized Case 

Recently, it was noticed that the BGF events in non-polarized DIS can be studied without 
involving the jet algorithms || . For this one can measure a liner interdependence between 
the current- and target-region multiplicities in the Breit frame. The approach involves 
the measurement of the particle multiplicities in large phase-space regions, rather than 
clustering separate particles at small resolution scales. Therefore, high-order QCD and 
hadronization effects are expected to be smaller than for the dijet studies that use clus- 
tering algorithms with specific resolution scales. In addition, the method does not depend 
on the jet transverse energy Et used in jet reconstruction and is well suited for low Q 2 
regions where the jet algorithms are less reliable. 

2 As usually, Q 2 — —q 2 = — (ft — k ) 2 (k and k denote the 4-momenta of the initial and final-state lepton, 
respectively) and the Bjorken scaling variable x is defined as Q 2 /(2P ■ q), where P is the 4-momentum of 
the proton. 
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The correlation between the current and target region multiplicities can be measured 
with the covariance 

cov = (n c n t ) - (n c )(n t ), (1) 

where n c (rat) is the number of final-state particles in the current (target) region and (. . .) 
is the averaging over all events. If h particles from the hard QCD processes are emitted 
in the target region, one can rewrite (|l|) as 

cov = ((n c - h) (fit + h)) - (h c - h)(n t + h). (2) 

Here h c is the total number of particles coming from zero and first-order QCD processes 
(h < ra c ) and fix is the multiplicity of the proton remnants without counting the particles 
due to the hard QCD processes. From (0) one obtains 

cov = (h c h) - (h c )(h) - (h 2 ) + (hf. (3) 

In this expression, the contribution from the remnant multiplicity ra t cancels since we 
consider the case when fit is independent of ra c and h. This key assumption means that 
the only important effects leading to the correlation between ra c and ra t are the hard QCD 
radiation. The validity of this assumption has been tested in Q for a Monte Carlo model 
based on the parton shower, describing higher than first order QCD effects, and the LUND 
string model for hadronization. 

At small x, the BGF is the only dominant first-order QCD process. Let us define the 
BGF rate as 

i?BGF = ^, (4) 

where -/Vbgf is the number of the BGF events and iV ev is the total number of events in 
the limit N ev — > oo. According to the assumption discussed above, h = for the QPM. 
For the BGF, however, one or two quarks are emitted in the target hemisphere, so that 
their parton radiation gives h > 0. The averaging in (|l[j^) is performed over all possible 
DIS events, despite the fact that cov ^ for the BGF events only. Therefore, one can 
apply the averaging over the relevant BGF events using the relation: 

(f(h)) = Rbgf (f(h)) B GF, (5) 

where f(h) is any function of h so that f(h = 0) = and (. . .)bgf is the average over the 
BGF events. Using this relation, one has 

cov = ((n c /l) B GF - ("c)(^)bGF - {h 2 )BGF + -RbGfWbGf) ^BGF- (6) 

The last term can be neglected since it contains a small contribution. Therefore, in the 
linear approximation, the covariance can be expressed as 

cov ~ —A(Q 2 ,x) Rbgf (7) 

with 

A(Q 2 , x) = (n c ) (/i) BGF + (/i 2 )bgf - (n c h) BG F- (8) 

A more detailed form of this expression has been obtained in ||. 
The function A(Q 2 ,x) depends on: 
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1) A number of the final-state hadrons in the jets initiated by the quarks in the BGF 
processes. 

2) Kinematics of the BGF jets in the Breit frame, which depend on Q 2 and x. Using 
LEPTO Monte Carlo @ with the tuning described in and GRV94 HO Q parame- 



terization of the parton distribution functions, we have found that for 5 < Q 2 < 50 GeV 2 , 
the fraction of BGF events with both quarks moving to the target region increases from 
52% at (x) ~ 0.2 • 10" 1 to 61% at (x) = 0.5 • 10" 3 (100% corresponds to all possible jet 
configurations of the BGF events in the Breit frame). Therefore, A(Q 2 , x) ~ A(Q 2 ) in (Q) 
is a good approximation. Note an x-dependence of the cov due to the BGF kinematics is 
rather small compared to the x-dependence of the BGF rate itself: -Rbgf increases from 
7% at (x) ~ 0.2 • 10" 1 to 21% at (x) = 0.5 • 10~ 3 . 

Before going to a polarized case, let us again remind all approximations made in (]?]). 

1. We neglect the QCD Compton scattering, considering low Q 2 regions. Note that even 
if a small fraction of the QCDC is present, some QCDC events cannot contribute to 
the correlations since singularities in the QCDC cross section favor the event topology 
with two jets in the current region, which does not produce the correlations (see, for 
example, Ref. 



2. Effects from high-order QCD and hadronization are already assumed to contribute 
to a value of A(Q 2 ,x), an exact form of which is beyond the scope of the present 
study. However, the assumption made to derive (0) was that the high-order QCD 
processes and hadronization do not change significantly the LO dijet kinematics in 
the Breit frame. This assumption leads to the factorization of A(Q 2 ,x) from the 
BGF rate in the linear approximation. 

Using a numerical simulation |J, it was shown that the effects quoted above do not 
strongly contribute to the relation (Jfy with an x-independent A(Q 2 ) in the range 5 < Q 2 < 
50 GeV 2 and 0.5 • 10 -3 < x < 0.2 • 10 _1 . According to the LUND model implemented in 
JETSET, the hadronization does not produce large correlations during the formation and 
independent breaking of the strings stretched between the current-region partons and the 
remnant. 



3 Polarized Case 

The ratio of spin dependent gluon density AG(x) to spin averaged gluon density G(x) at 
a fixed Q 2 is proportional to the asymmetry of the BGF cross sections, 

A±U{X) _ ObGF "BGF _ °tot -"-BGF "tot -"-BGF /q\ 
{ ' ^BGF + BGF a tot K BGY + a tot K BGF 



where (kill) is the value of the BGF asymmetry at partonic level [16|, cr^ GF (R^ GF ) and 
°bgf (-^bgf) are * ne BGF cross sections (BGF rates) in case of the ant-parallel (f J,) and 
parallel (ft) polarizations of the incoming lepton and nucleon. a^ t and o tot are the total 
cross sections which can be obtained by counting the number of events in a given kinematic 
bin for the different spin configurations, normalized to an integrated luminosity. 
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The expression (0) can further be rewritten as 



(a^^^^Ap + Au, Ap = ^ , (10) 

BGF + BGF 



n _ tt 

A °tot °tot / 1 1 \ 

°tot + "tot 

where An is an inclusive polarized asymmetry, which is expected to be small compared to 
the asymmetry Ap of the BGF rates. The relationship between @ and ( |l0|) can easily be 
seen after rewriting Ap + An as 

°tot ^BGF °tat -"-BGF _ _ V^BCF n BGFA°toi ° tot) 



l^tot -"BGF + °tot -"bGFA 1 - a ) l \ a tot ^BGF + °tot ^BGfJ 

which is approximately equal to (||) by neglecting a small contribution from a. 
According to (|7j), 

cov U „ _ a( q2) i? U QF! CQV TT „ _ A( q2 } ^TT gf! (12) 

where A(Q 2 ) is considered to be independent of the polarization, since this function is 
mainly determined by the number of the final-state hadrons in jets, i.e. by multiple-gluon 
radiation and hadronization. From (12) it is easy to see that the theoretical Ap can be 
determined via the covariances 

. cov^ — covTT 

Ap ~ ACOV = T . rqr. (13) 

covTi + covtT v ; 

Thus the asymmetry AG(x) /G(x) can experimentally be obtained from Acov and directly 
measurable A\\. Below we shall numerically estimate the BGF asymmetry Ap (|To|). We 
calculate this quantity directly from the BGF rates and then reconstruct it by measuring 
the Acov from the final-state hadrons of polarized DIS. 



4 Numerical Studies 

In order to study the asymmetry using the current-target correlations we use the PEPSI 



Monte Carlo generator [17] for the polarized leptoproduction. This model is based on the 
LEPTO 6.5 for DIS together with JETSET 7.4 describing the fragmentation. Two DIS 
samples with the opposite spin configurations were generated in the region 5 < Q 2 < 50 
GeV 2 . The electron and proton beam energies were taken to be 27.5 and 920 GeV, 
respectively. For simplicity, both beam polarizations were assumed to be 100%. The total 
number of the generated DIS events is 1M for each polarization sample for the given Q 2 
range. The covariances for each polarization were determined from charged final-state 
hadrons according to (||). Hadrons with lifetime cr > 1 cm are declared to be stable. 

Fig. |2|a shows the behavior of Acov calculated from final-state charged hadrons com- 
pared to the theoretical prediction for Ap, generated also with the PEPSI using informa- 
tion on the generated type of the LO process. The GS-A f\M polarized parton distribution 
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Figure 2: The asymmetry calculated using correlations between current- and target-region 
charged multiplicities generated with PEPSI (full symbols) compared to the theoretical ex- 
pectations for Ap of the same model (solid lines). The figures show the x-dependence of 
the asymmetry in the range 5 < Q 2 < 50 GeV 2 : (a) The asymmetry for the GS-A model 
implemented in PEPSI; (b) The same but when AG is set to zero. 



was used. From this figure one can conclude that the experimentally measured asymme- 
try Acov does reproduce the theoretically expected behavior of the Ap. This is especially 
evident for a small x, where the contribution from the QCDC can be neglected. It is also 
seen that typical statistics expected at HERA are sufficient to reliably determine the sign 
and the magnitude of the asymmetry. 

For an illustration, Fig. shows the same as Fig. ^a, but after switching the gluon 
polarization term off, i.e. AG = 0. The quark polarization distributions were unchanged. 
In this case, both Ap and Acov tend towards zero. A difference between Ap and zero 



value is expected to come from a positive value of .Am (11). It may also be noted a small 
discrepancy between Acov and Ap. This might mean that for this, in fact, unrealistic 

1" I TT 

case, the values of R^ GF and R^ GF are so close to each other that an approximate nature 
of the expressions fll2] ) makes a noticeable effect on the relationship between Acov and 
Ap. 

Dependence on cut-offs. The LEPTO Monte Carlo used by PEPSI contains cut-offs 
to prevent divergences in the LO matrix elements. Therefore, the magnitude of the BGF 
rate in PEPSI is ambiguously defined. This may produce a systematic bias in the relation 

Let us remind that, in the BGF cross section, the singularities in the two-parton 
emission are given by l/z(l — z) |0J, where z = (jp' -p)/(p' ■ q) (p is the momentum of the 
incoming parton and p 1 is that of the final-state parton, q 2 = —Q 2 ). The LEPTO cut-offs 
are based on the so-called "mixed scheme" in which the parameter z m i n , restricting the 
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values of the variable z, plays an important role in the determination of the probability for 
the BGF event to occur. The default value of this cut-off is set to 0.04. If one decreases this 
cut-off, the relative contribution of the BGF rates increases, with respect to the QCDC. In 
this case, the relation ([?]) is expected to be a good estimate. However, for a large z m i n , the 
BGF rate is small and Acov might poorly reflect the behavior of Ap. Therefore, for our 
systematic checks, the value of z m i n was increased. We have verified that the relationship 
between Acov and Ap still holds up to z m i n = 0.2 for the given statistic uncertainties. 
Note that while for the default value z m i n = 0.04 the BGF rate is about 0.6 in the smallest 
x-bin, the rate of this process drops down to 0.3 when z m i n = 0.2 is considered. 

We have also found that Ap changes only a little by varying the cut-off. This can 
be explained by the normalisation (-Rbgf + -^bgf) use d in (0)- Note that, from the 
experimental point of view, a similar normalisation (cov^ + cov^T) in would help to 
measure the Acov reliably even if the detector track acceptance for the target-region is 
small (see [11| for details). 

Dependence on the structure function. An important test for our method is to inves- 
tigate the asymmetry measured with Acov for different types of the parameterization for 
the parton densities. Various types of polarized parton distributions included into the 
PEPSI package were analyzed. The method was found to be sensitive to input structure 
functions and Acov reproduces the trends of Ap. 



5 Conclusions 

In this paper we have proposed a new method to measure the asymmetry in polarized 
lepton-nucleon scattering. This method is based on the measurement of the current- 
target correlations in the Breit frame. The advantage of the method is that all inclusive 
DIS events with well reconstructed Breit frame can be used to determine the asymmetry, 
without specific constraints to select useful events for this measurement. 

In this respect, it is important to emphasize that our approach is well suited for rather 
low Q 2 and Et, i.e. for the regions where dijet reconstruction suffers from misclusterings 
and large hadronization corrections. Thus the suggested method compliments and extends 
the study of AG/G to kinematic regions of low Q 2 (Et) where the method suggested in @ 
is less reliable. We also expect different systematics for these two methods: While the 
method discussed in is a subject of the systematic uncertainties as for the standard 
dijet studies, systematic effects for our method most probably would come from a low 
detector acceptance in the target region of the Breit frame. 
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